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Abstract 

We demonstrate how one can describe explicitly the present arbitrariness in 
solutions of relativistic wave equations in external electromagnetic fields of 
special form. This arbitrariness is connected to the existence of a transforma- 
tion, which reduces effectively the number of variables in the initial equations. 
Then we use the corresponding representations to construct new sets of ex- 
act solutions, which may have a physical interest. Namely, we present new 
sets of stationary and nonstationary solutions in magnetic field and in some 
superpositions of electric and magnetic fields. 

I. INTRODUCTION 

Relativistic wave equations (Dirac and Klein-Gordon) provide a basis for relativistic 
quantum mechanics and quantum electrodynamics of spinor and scalar particles In 
relativistic quantum mechanics, solutions of relativistic wave equations are referred to as one- 
particle wave functions of fermions and bosons in external electromagnetic fields. In quantum 
electrodynamics, such solutions allow the development of the perturbation expansion known 
as the Furry picture, which incorporates the interaction with the external field exactly, 
while treating the interaction with the quantized electromagnetic field perturbatively 0. 
The physically most important exact solutions of the Klein-Gordon and the Dirac equations 
are: an electron in a Coulomb field, a uniform magnetic field, the field of a plane wave, the 
field of a magnetic monopole, the field of a plane wave combined with a uniform magnetic 
and electric fields parallel to the direction of wave propagation, crossed fields, and some 
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simple one- dimensional electric fields (for a complete review of solutions of relativistic wave 
equations see [Q). 

Considering, for example, stationary solutions of relativistic wave equations, we can see 
that in the general case, there exist different sets of stationary solutions for one and the 
same Hamiltonian. The possibility to get different sets of stationary states reflects the ex- 
istence of an arbitrariness in the solutions of the eigenvalue problem for a Hamiltonian. 
Considering nonstationary solutions, we also encounter the possibility of constructing dif- 
ferent complete sets of such solutions. There is no regular method of describing such an 
arbitrariness explicitly. Especially in the presence of an external field the problem appears 
to be nontrivial. 

In the present article we demonstrate how one can describe explicitly the present arbi- 
trariness in solutions of the relativistic wave equations for some types of external electro- 
magnetic fields, namely, for uniform magnetic fields and combination of these fields with 
some electric fields. This arbitrariness is connected to the existence of a transformation, 
which reduces effectively the number of variables in the initial equations. Then we use the 
corresponding representations to construct new sets of exact solutions, which may have a 
physical interest. In Sect. II we consider relativistic wave equations in pure uniform magnetic 
fields. Here we derive a representation for the exact solutions, in which the above mentioned 
arbitrariness is described explicitly by an arbitrary function. From a suitable choice of this 
function, we get both the well-known set of solutions and new ones. This Section contains 
the most complete (at the present) description of the problem of a uniform magnetic field in 
relativistic quantum mechanics. Among new sets of solutions there are both stationary, gen- 
eralized coherent solutions and nonstationary solutions. Then, in Sect. HI, we consider more 
complicated configurations of external electromagnetic fields, namely, longitudinal electro- 
magnetic fields. Here we describe all the arbitrariness in the solutions, and on this base 
present various sets of new exact solutions. In Sect. IV we interpret the above results from 
the point of view of the general theory of differential equations. 



II. UNIFORM MAGNETIC FIELD 

A. Arbitrariness in solutions of relativistic wave equations. 

Consider a uniform magnetic field H = (0, 0, H) directed along the axis {H > 0). The 
electromagnetic potentials are chosen in the symmetric gauge 

^0 = ^3 = 0, Ai = ^Hx^ A2 = ~Hx'. (2.1) 
We write the Klein-Gordon and the Dirac equations in the form 

/C* = 0, fi^lC = p2 _ Vf, = ihd^ - -A^ , 

= , hV = 7^P^ - moc . (2.2) 

Here e = — |e| and 7-matrices are chosen in the standard representation [Q. 

In the field under consideration, the operators Vq and are mutually commuting inte- 
grals of motion, [/C,Po] = [/C.Pa] = [V,V^] = \V,V^] = [Vo^V-^] = 0. 
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In the case of the Klein-Gordon equation, the operator L^, 



(2.3) 



can be included (together with Vo and V3) in the complete set of integrals of motion, whereas 
for the Dirac equation case, the operator J^, 



J, = + ^Sa, [J„ Vo] = [J„ Vs] = [V, J,] = , 



(2.4) 



can be included (together with Vo and Vz) in the complete set of integrals of motion. Here 

E3 = diag (0-3, as) . 

We are going to use dimensionless coordinates —00 < x < 00, —00 < y < 00 or 
0<p<oo, 0<(/3<27r defined by the relations 



^x'^ = x = ^cos^, ^x"^ = y = y/psimp , 7 = > , 
2 1 

dx^dx^ — —dxdy — —dp dip, x + iy — A/pexpiip . 
7 7 

It is useful to introduce the operators ai, 02, Oi", 0^2^ , 



(2.5) 



02 = - 


1 r 


s/2^h \- 


«2 = - 


1 r 




Oi = - 


1 




4 = - 


' ( 





V2 - iVi + h'j (x^ + ix^) 



1 



^.x + iy + d^ + idy) - 



{p + id^ + 2pdp) , 



—lip 



V2 + iVi + h'j (x^ - ix"^) = -{x-iy-d^ + idy) = — — (p + id^ - 2pdp) 
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[iVi +V2) = -{x-iy + d^- idy) = ^-^ (p - id^ + 2pdp) , 



-'P2) = ^{x + iy - 
They obey the commutation relations 



3«V 



2Vp 



(p - ia,^ - 2pdp) . 



0, A;,s = l,2. 



(2.6) 



(2.7) 



Thus, we can interpret these operators as creation and annihilation ones. One can also find 
the following relations 



(2.8) 



Pi + V2 = [aia+ + a+aij = 2h^'yAf + ^^7 , 
Lz = h (Af — a2 02) , Af = afai . 

Then the Klein-Gordon and the Dirac operators can be written as 

IC = h-'(v',-Vi)-2jM-j-m', m = 
V = h-' (7°7'o + l^Vs) - [(7' - ^7') «i + (7' + ^7') at - 



m . 



(2.9) 



The operator TV commutes with Vo, V^, Lz, plus it is an integral of motion in the 
case of the Klein-Gordon equation. Its generalization for the Dirac equation has the form 

One ought to remark that the operators /C and V do not contain the operators 02- 
Thus, the latter operators are integrals of motion, which commute with J\f, Afn, Vq, V3, but 
do not commute with and Jz ■ 

The operators of creation and annihilation with different numbers commute. One can 
find a representation in which these operators are acting on different variables. To this end, 
we present the wave functions from ( |2.2| ) in the following form (we make a Fourier transform 
in the variable y only, and call such a representation the semi-momentum representation) 

^(x, y) = ^= / e^'^y^^ix, k) dk . (2.10) 
V 27r "'-00 

Of course the functions \1/ and \1/ depend on the variables x'^ and x^ as well, but we do 
not indicate this dependence explicitly. In terms of \I' the multiplication and differentiation 
have the form y —>■ idk, idy — —k. Then, the expressions for the creation and annihilation 
operators in the semi-momentum representation take the form 

2ai = X + k + dx + dk, 2a^ = x + k — — dk , 

2a2 = X — k + dx — dk, 2a^ = x — k — dx + dk ■ (2-11) 
Now we pass from x, k to new variables ^, 77, 

V2^ = x + k, V2r] = x-k, V2x = ^ + r], V2k = ^ - r] . (2.12) 
Then the creation and annihilation operators can be written as 

V2ai = ^ + d^, V2at = ^-d^, V2a2 = v + drj, \f2at = i]-d^. (2.13) 
In the new variables, 

2M = -dl-\, (2.14) 
and the Klein-Gordon and the Dirac operators read 

}c = n-' {vl-vt)^i(dl-e)-^\ 

V = h-' (7°Po + I'Vs) - ^ (fi - i^^d^ - m . (2.15) 

One can see that the latter operators do not contain the variable rj. Notice that both 
operators and Jz contain variables S,,ti. For example, 

2Lz = e-dl-r,^ + dl. (2.16) 

The integration over k in ( p^.lO[ ) can be replaced by an integration over r], 

00 

^{x,y) = ^ e~'-^y'^^{^,v) d7], ^ = V2x - v ■ (2.17) 

—00 
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Besides, one can write 

oo oo 

(^,$)= J dx J dy^*{x,y)<^{x,y) = {^,^)= J J dv ^* {^,v) ^ {^,v) ■ (2.18) 



oo oo 



-oo — oo 



-oo — oo 



The independence of the operators ( [2.15| ) on the variable t] will allow us to separate 
explicitly the functional arbitrariness in the solutions ( p.l7| ), as will be seen below. 



B. Stationary states 



Known sets of stationary solutions in a uniform magnetic field (that were found in the 
first works are eigenfunctions of the operators Vo,V3,Af in the scalar case and of 

the operators Vo, V3, Md in the spinor case. Thus for scalar wave functions ^ we have the 
conditions 



1 



n = 0,1,2, 



Vq^ = hko"^, = hks"^ , A/"^ = n^, n = 0, 1, 2 

and for Dirac wave functions ^ the conditions 

Consider first the scalar case. It follows from ( p. 15 ) that 

fcg = fn'^ + 7 + ^3 + = m* ^ + /Eg + 2^n, m* ^ = + 7 

and 



,,fc3 {x^") = Nexp [-ikox^ - ik^x^) ^„ (x, y) 



(2.19) 



(2.20) 



(2.21) 



(2.22) 



Here is a normalization factor. In the semi-momentum representation ( p.lO|) the function 
^ra (x, y) has the following image 



X — 7] 



(2.23) 



Here Eqs. (|2.19| ), ( p.l4|) were used, f/n (0 are Hermit functions; they are related to the 

corresponding polynomials -ffn (0 as Unix) = (2"n!y7r)~^ exp (— x^/2) if„ (x) 
function $ [r]) is arbitrary. The functions [x, y) from (|2.22|) obey the relations 



ai^„ = v^^„_i, a+^„ = Vn + l^n+i , (x, y) 



r(n + l) 



00 

(^5 y) = exp + ixy^ j drj exp 



-— + V2r7(x - iy) 



The 

(2.24) 
(2.25) 



Dirac wave functions are of the form "^n^i^^ 
with bispinors ^ks,n {x-, y) having the structure 



$(r/) . 

exp {-ik^x^ - ik^x^) ^^g,™ {x, y) 
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*I,fc3(^'?/) = (ci^n-i(a;,i/), ic2'^n{x,y), C3^'„_i(x, ?/), ic4'i!nix,y)) . 



(2.26) 



The functions \I/„(x,?/) are defined by the relations ( p.l7|) , ( p.23|) , whereas the constant 
bispinor C (with the elements c^) obeys an algebraic system of equations 



The condition det A 
of {^21 



AC = 0, A = 7°A;o + -f^h - ^^1^ - m . (2.27) 
= (/cq — kl — 2'~fn — m?Y — results in an equation which is an analog 



ftp 



^3 + 27n + 



f2.28) 



Since the rank of the matrix A is equal to 2, a general solution of ( |2.27| ) has the form 



C 



(/cq + m)v 



C+C = 2kQ{kQ + m)v^v 



(2.29) 



where v is an arbitrary constant bispinor and cr are Pauli matrices. We can specify v 
selecting a spin integral of motion (see 0). The state n = is a special case. Here we 
must set Ci = C3 = 0, that corresponds to the choice v'^ = (0, C2) , C2 7^ 0. The latter 
means that Ti^'^d = —"^d- Thus, for n = 0, the electron spin can only point to the direction 
opposite to the magnetic field. 

Expressions for \E'„(a;, y) in the semi-momentum representation contain explicitly a func- 
tional arbitrariness, which means that every energy level is infinitely degenerated. Let us 
demand that the scalar and spinor wave functions be eigenvectors of the operators and 
Jz respectively. According to ( p.4|) and ( p.8|) that means that the functions \E'„(x, y) have to 
obey an additional condition 



ata2'^n{x,y) = s'^n{x,y), s = 
= h (n — s) = hL I = n — s. 



0,1,2,... 



n > I > —00, X 



(2.30) 



This condition defines the function $(27) according to ( p.l3|) , a2a2^s{v) = ^^s{v)y therefore 
$s(^) = Usiji). Substituting this result into ( |2.23| ) and into ( |2.17| ), and doing the integral 
over 7], we find in the coordinate representation. 



^n,s(a;,i/) 



[—1)^ f X + iy^ 



'2tt \x-iy 



(2.31) 



Here Im,n (x) are Laguerre functions, which are connected to the corresponding polynomials 
(x) by the relations (see [|1^) Im,n (x) = (F (n + 1) /T (m + 1))^^^ e~^x^L" (x) , a = 
m — n. The states ( p. 31 ) were first obtained in the works Besides ( p. 24 ) and ( p. 25 ), 

the functions (|2.31|) obey the following relations as well 



a2^ 



Vs + 1* 



n,s+l ; 



n,s 



r{n + i)r{s + i) 



^0,0 , ^o,o(a;,l/) = — exp 

TT 



e 2 



(2.32) 
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Below we are going to find new sets of solutions imposing complementary conditions 
different from ( 2.30|) . This results in a different form for the function ^{r]). 

Taking into account that the operators 02", 02 are integrals of motion, we may construct 
stationary states, which are eigenvectors of a linear combination ^2'^ of these operators. 



(2.33) 



Here are arbitrary complex numbers. One has to distinguish here three nonequivalent 
cases: 

If lap < then do not exist any normalizable eigenvectors of the operator (|2.33 ) . 
We are not going to consider such case. 

If |ap = then A"^'^ is, in fact, reduced to a Hermitian operator 



Ai^ = ^,a2 + ^i*at, At^ = A^2, /x ^ 



(2.34) 



where /x is an arbitrary complex number. 

If lap > |/3p, then without loss of generality we can assume that operators ^2'^ have 
the form 



^2'^ = aa2 + (^0,2^ 



a 



^2 ) ^2 



1. 



(2.35) 



Then are creation and annihilation operators, which are related to 02", 02 by a 

canonical transformation 



02 = aM^'^ - pA. 



1+ a,l3 



aAf''^ - f3*A^''' . 



(2.36) 



Consider eigenvectors of the operator ( |2.34]) , i.e., Ai^'^i^ ^{x,y) = ^(x, y), z = z*. This 
equation results in the equation ^2$^ (77) = 2;$^ (77) for the function ^{rj). Taking into 
account ( p.l3| ), one can find that solutions of the latter equation are 



fx 



expQi 



V^TT |/i| (/i — /i* 

4{fi- fi*) Qi = -2 (/i + fi*) + AV2zrf - (/i + /i*) l/if' . 
These solutions obey the orthonormality and completeness relations 

/OO POO 
^yiriWAv) dv = 6{z -z'), / $: ''WWM dz = 6iv - V') 
-00 J —00 

Their overlapping has the form 



(2.37) 



(2.38) 



(^'^ z) = 1 <^*/'{vM ^iv) dv = Ni exp 

Q2 



Q2 



4 (/i'/i* — nn'*) 



/i /i 

27?^ l/i'l |/i| (/i/i'* — /i'/i*) 



/i 




(2.39) 



It defines the mutual decomposition 

oo 

$^ [r]) = J [t]) R^''^ {z', z) dz' . (2.40) 

— oo 

The coordinate representation ( p. 17 ) for the solutions under consideration has the form 



* \ 2 



vl/^:^, {x,y) = (v^TT l^^j ' Un {pi)exptQ, 

4 |Ai|^ Q3 = [i (/i* - fi)x + {fi + fi*) y] [{fx + n*)x + fi*)y - 2z] , 

\/2 = (/U + /U*) X + i{fi — fi*) y — z . (2-41) 

Their scalar product ( |2.18| ) reads (^^n',z'^^n,z) = ^n,n'S{z — z'). The relation ( |2.4CI| ) results 
into the following decomposition in the coordinate representation 

00 

K,. y)= j Kz' y) R^'^' (^'' ^) dz' . (2.42) 

—00 

In particular, in the cases of real or pure imaginary /x, such wave functions were known 
before 0. 

Consider eigenvectors of the operator (|2.35| ), i.e., A'^'^'^^z i^yV) = ^"^nz {^^v) where z 
is a complex number. In fact, we get coherent (squeezed) stationary states. They are labeled 
by z and by two complex parameters a, (3, which are related by the condition (|2.35|) . In the 
semi-momentum representation the above equation is reduced to the one 

{V) = iv) ■ (2.43) 

It is well known that such solutions form a complete (overcomplete) set at any fixed a, j3. 
Solutions within each set are not orthogonal. One can use these functions to construct a 
orthogonal set of solutions. 

Since the operators A2 "^'^ , are integrals of motion (both for the Klein-Gordon 
equation and for the Dirac equation), they are symmetry operators for the equations. The 
action of these operators on a solution provides again a solution. For example, applying 
the operators (F (1 + s))'^^"^ (^A^ "'^ — z*^ , s = 0, 1, 2, ... to normalized solutions of the 
equation ( p.43|) , we get normalized solutions labeled by the index s. These new solutions are 
orthogonal with respect to s, 



T{l + s) 

00 

/ (r/) $:;f iv) dv = Ss,s' J \'^:'^ iv)\' dr^ . (2.44) 

—00 —00 

We call such states generalized squeezed coherent states. It is possible to get an explicit 
form for these states, 
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a 



la I (a — f3) 



'a* - (3* 



2 {a/3* - a*P) + 2^27] [z {a* - /?*) - z* {a - + z*^{a- 

-z^ (a* - p*f , 2 |a - /?| p2 = 2r/ - (a* - /?*) - V2z* («-/?) • (2.45) 



(2.46) 



The functions (|2.45|) form a complete set for each fixed z , 

00 

E*:;'"(V)$:;f(^) = 5(v-r/), 

and for each fixed s , 

72 

/ — <l>:"'^(r/')Cf(^) = <5(^-^?'), rf2^ = rfRezdIm2. 
The overlapping 

J —00 

allows us to find mutual decompositions 

00 „ 
^tUv) = E <f •"''(^', ^Tfiv) = / d'z' Rif'-'^iz', ^)$^;i!'(r^) . (2.49) 

s'=0 



(2.47) 



(2.48) 



Unfortunately, the overlapping (|2.48|) has a complicated form via a finite sum of Hermit 
functions. In some particular cases this sum can be simplified. For example, if a' = a, 13' = 
P, then the overlapping does not depend on a,P and has the form 



i?:;J"'^(z',^)=i?,,,(z',z 
For s = s' = we get 



Z — Z 



exp 



,zz — z z] 



/,,,(|^-zf). (2.50) 



■0,0 



-{aa'* - 1313'*)'^' expQ 



5 ; 



2Q5 



z^(a^*/3* - a*P'*) + (z^*)^(a/3^ - a^/3) + 2zz'* 
aa'* - (313'* 



z z — zz — \z — z 



'* |2 



(2.51) 



The wave function \E'"'f ^(x, y) has also a very complicated form in the general case. However, 
in some particular cases it can be simplified. For example. 



Vvr \x — ty — z* 
M = z{x — iy) — z*{x + iy) . 



(2.52) 



For z = we arrive at the set (|2.31|) . 

For s = we get a compact form for a set of stationary squeezed coherent states 
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P3 



exp Qe 



\a\ \^) " \ V^2a/?^ 
p^ = z-a{x + iy) - - iy) , 

AapQe = (1 + 2|/3|2)z2 _ 2a(3\z\'' + {z + ps) [p{x - ly) - a{x + iy)] 
Additional simplifications are available for a = 1, /5 = 0, 



(2.53) 



(x + - 2)" 
^nA^^y) = I ^exp 
Vrfn + 1) 



z(x -iy) - - (a;^ + y^) 



(2.54) 



Namely, these states where found in 0. However, the meaning of the parameter z was not 
clarified. 

For arbitrary a, /?, the functions '^'^''^ ^{x,y) obey, besides (|2.24|), the following relations 



(2.55) 



Taking into account that a2ipn,z = dipn^z/dz, we can construct a new set of stationary states 
by successive differentiations. 



^n,s,z ("^j y) 



'IT 



■ exp 



exp 



i{n — s)ip + 



{x + iy) 



p-q 



1 



n — s 
2 



X + iy — z 
X — iy 



n — s 
2 



h,n{q) 



q = p — Zy^e ^'^ = {x — iy){x + iy — z) . 
For = 1 the above set obeys (besides (|2.24 )) the relations 

d 



02^ 



n,s,z ^ ^ n,s,z 



Z'^r 



y/s'^n,s-l,z , at'^n,s,z 



dz 



(2.56) 



(2.57) 



The set ( |2.31D is a particular case of ( |2.56D , it corresponds to z = 0. The set ( |2.5(j| ) is not 
orthogonal. 



[^n',s',z','i'n,. 
Js,s' {Z, Z) = 

J.s' (Z, Z) = 



N'*N6n,n'Js,siz,z') 
Z^'-s^zz'"^ ' 



r 


(^ + 1) 


r( 


.' + 1) 


r( 


.' + 1) 



S'-S„ZZ-TS'-S(^^,*. 



s<s', 



\ r(. + i) 



(2.58) 



The functions from the set ( p.56| ) are normalized to unity for N = Ns{z) 
exp {—\z\'^/2) [Lg (— |zp)] For N = 1, the following mutual decompositions take place 
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n,s+k,z' 



T{k + 1) f£z 



\T{s + k + l) J TT 



n,k,z 



T{k + s + l){z' - z) 



k\ 



n,s+k,z 



(2.59) 



That means, in particular, that ( |2.56| ) is a complete set since the set ( |2.31| ) is complete. 

Selecting different forms for the function $ (77) , we can get other sets of stationary states 
for a charge in a uniform magnetic field. 



C. Nonstationary states 

The most interesting nonstationary solutions of relativistic wave equations for a charge in 
a uniform magnetic field are coherent states; for the first time such solutions were presented 
in [p!0|-p!3[], see also 0. Below we present a new family of nonstationary solutions, which 
includes the above coherent states as a particular case. 

Here we are going to use light-cone variables vP = x'^ — x^, = x^ + x^, and the 
corresponding momentum operators 

Vo = ^h^o = li^o- V3) , Vs = thds = lin + Vs) , (2.60) 

where do = d/dvP, 83 = d/du^. Then the Klein-Gordon operator can be presented in the 
form 

K. = Ah-^VsVo - -m*\ (2.61) 
whereas the Dirac equation reads (^ is a Dirac bispinor) 

4r2p3^o^(-) = (27A/'d + , 2P3^(+) = [iaV±) + kp^rn] ^(_) , 

V^ = -{Vi,V2,0) , * = ^(+) + ^(_), ^(±)=p±*, 2p± = l±a3. (2.62) 

Here a and pa are Dirac matrices [Q, and p± projection operators. 

In the case of the uniform magnetic field under consideration, the operators V^yVo are 
integrals of motion. Thus, we will consider solutions that are eigenvectors of V3, 

= . (2.63) 



The scalar wave function obeys ( |2.63|) and can be written as 



. A .3 . m 



* 2 

^[x^) = A^exp ( -i-u^ - I'-^u^ ] ^(u°, X, y) . (2.64) 
\ 2 2A / 

It is easy to see that ■ip{u'^,x,y) obeys a first order equation, which can be treated as a 
Schrodinger equation. 



idoip (^u^,x,y'^ = uafaiip (^u^jXjy^^ , u = — . (2.65) 
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Suppose Eq. (|2.63|) holds, then can be presented in the form: 

/A m*^ \ 

^(_)(a:'^) = N exp i-i-u^ - i^u^ \ W {1 - 03) C^{u^,x,y) . (2.66) 

Here C is an arbitrary constant bispinor, and is a unitary matrix [ipo is a constant phase), 

= cos K - iSg sin K, 2k = uju^ + (fo, W+W = I , (2.67) 

and ijjivP^x^y) is a scalar function. The latter function obeys the equation ( |2.65| ). Then, 
the projection can be found from ( f2.62D , ^{+) = {^^)~^ [('^'^±) + ^"^Ps] ^(-)- 

Thus, both in the scalar and spinor cases we have to solve the same equation ( p.65| ). 

In the semi-momentum representation, the corresponding function ip{vP,^,ri) obeys the 
same equation ( |2.65| ), where, however, one has to use the expression ( |2.14| ) for the operator 
M = afai. The relation between the functions ip{u^,C,,Ti) and ip{u^,^,r]) still has the form 

Let us introduce the operators 

A{'' = fa,+gat, ^'^ = fat + g*a, , (2.68) 

where the complex quantities / and g can depend on u^. These operators are integrals of 
motion whenever /, g obey the equations (by dots above are denoted derivatives with respect 
to M°) 

if + ujf = 0, ig-ujg = 0. (2.69) 

It is easy to find 

/ = /o exp {iuju°^ , g = go exp (^-iuju'^^ , (2.70) 

where /o, go are some complex constants. Bearing in mind considerations related to the 
operators (|2.33|) , we are going to consider two nonequivalent cases only. The first one 
corresponds to |/p = or equivalently to |/oP = IfifoP- In this case we can, in fact, only 
consider the Hermitian operator 



= uai + u*at, V = z/oe*'^"°, z/q = const . (2.71) 
The second case corresponds to |/p > \g\'^, and here we can suppose that 

\f\'-\g\' = \fo\'~\go\' = i, (2.72) 

without the loss of generality. In both cases the operators ( |2.68| ) are, within constant complex 
factors, creation and annihilation operators. 

Let us include operators (|2.71|) and ( p. 341) (they are integrals of motion) into the complete 
set of operators. Then 

A'^r.[% = ^ir.[%, Ai^r.[% = Z2r.[%, 4 = ^k, k = i,2. (2.73) 

In the semi-momentum representation we find 
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^{u',^,r^) = ^:^{OKiv) , (2.74) 
where functions are defined in ( p.37| ). The corresponding coordinate representation reads 



exp 



4 (/i*z/* — /iz/) 



27r2|/i||z/|(yuz/ — /i*z/*) 
Qe = 2(/i + ^J*){v + z/*)x^ + 2(/i - /i*)(z/ - z/*)?/^ + 4i(/iz/* - /iV)x?/ 
-4x [2:1 (/i + /i*) + Z2(z^ + T^*)] - ^iy [ziifi - /i*) - ^2(z^ - z^*)] 

rjjF fir \ ( nr ^ 



(2,75) 



These solutions are orthogonal at any fixed m", 
obey the completeness relation 



=S{z,- z[) 5 {z, - 4), and 



dz, / dz2 r/JXu\ x', y'W.t.Sn\ x,y)=b{x- x') S (y - y') . 



(2.76) 



Consider now generalized squeezed coherent states, which can be constructed by anal- 
ogy with (|2.44|) in the semi-momentum representation. We use here the operators (|2.68|) 
supposing that the relations ( p39|) , (pTO| ), (^), and (^) hold, 



The functions (x) are defined in ( ^.45] ). Thus, 



(2.77) 



^2 ^2 



r(n + i)r(s + i) 



(2.78) 



The solutions ( [^.78| ) obey the relations 



21,22 ' 



(2.79) 



Eq. ( p.78|) describes the most general form of relativistic wave equation solutions in a 
constant uniform magnetic field. All the formerly known solutions can be obtained from 
this equation by a particular choice of parameters. For instance, by selecting /o = « = 
1, g = [3 = zi = 0, Z2 = z with 2; = we get the states ( p^.31D , on the other hand, if one 
puts s = 0, z 0, then one gets the states ( p.54|) . For n = s = 0, fo = a = 1, g = P = 0, 
we get coherent states [T^-|T3[. 

In the general case, an explicit coordinate representation for the solutions (|2.78|) looks 
complicated enough. However, some particular cases admit essential simplifications. For 
example, suppose /o = a = l, g = P = 0, then 
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1,0:1,0 



-1)"- fx + iy-zl-Z2^ 



X — ly — zi 



2Mi = {zi — Z2) {x + iy) — {zl — Z2) (x — iy) + z\z2 — ziz*2 — 2inuju^ , 

P4, = \x + iy — zl — Z2\'^ , zi = zi exp(— zdju") , (2.80) 

For n = s = 0, we get the coordinate representation for the squeezed coherent states in the 
form 



Q7= ^ 



af 



[af - l3g)'K\a 



zi\^ + \z2n + 



2 V'-- ■ y ■ 2iaf-/3g) ' 
q = -{a + (3)U + g)x' -(«-/?)(/- g)y^ + 2z(/3/ - ag)xy 
+2x [(« + P)z^ + (/ + g)z2] + 2iy [{a - - (/ - g)z2] 



Hag* - f^Dzf + if3*f - a*g)zl - 2z,Z2 . 



Solutions from [p!0| -[13| are particular cases of ( p.81|) for /o 
Calculating mean values in the states ( |2.78D , we getQ 



a 



1, g = P = 0. 



Vi = thJM{r + g*)z,-U + 



"ij ' 



V2 



(2.81) 



(2.82) 



Here we have taken into account the relations (|2.6| ), ( p.36|) , (|2.79|) , and the orthogonality 
of the states with respect to the indices n, s. Remember now that in classical theory the 
corresponding momenta , Vf have the following parametric representation (with vP being 
the evolution parameter, R radius of the classical orbit, and n is given by (|2.67|) ) 

= h-iR sin 2k, = -filR cos 2k . (2.83) 

It is easy to see that (|2]8|) coincides with jWM) for zi = (7/2)^/^ R {foe-'^° + goc'^^). Cal- 
culating mean values of the coordinates x^, x^, we find that they evolve as the corresponding 
classical quantities x^ 



,1 cl ^2 cl 



•^(0) ' ^(0) coordinates of the orbit center) 



x^^' = i^cosK + a;(o), x^"^ 



(2.84) 



for Z2 = (7/2)'/' [{a + (3)x\^^ + i{a - (3)x%^^ 

Thus, mean- value trajectories in the plane do not depend on quantum numbers 

n, s. These trajectories have classical forms under a proper choice of zi,Z2 . 

Calculating quadratic fluctuations in the states ( ^^.78] ), we get 



mVi) = ^'l\f + gVC^n + 1), 2(AP2)' = h'i\f-g\' (2n + 1) 
27(Aa;i)' = \f - g\\2n + 1) + |a - /3p(2s + 1) , 



27(Aa; 
o-i 



2Y = \f + g\\2n + l) + \a + (3\\2s + l) 
cj2 = tifg*~gni2n + l), 



ak = {Ax'^) (APfc) + (APfc) (Ax*^), A; = 1, 2 . 



(2.85) 



^One can obtain the same results using spinor wave functions for the calculations. 
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They do not depend on 2:1,-22, but do depend on quantum numbers n, s and on parameters 
fo,go,a,p. The relations ( p.85|) imply the generalized Heisenberg inequalities 



4Ji 



(2n 
{2n 



(2s + l)|(«-/?)(/ + ^?)p 
{2s + l)\{a + m -9)\' 
k = l,2 . 



(2.86) 



One can fix any given {Ax'') or {AVk) in a given "instant" by means of a choice of 
parameters fo,go,Ci,P ■ Then they evolve with "time" m° according Eqs. ( p.85| ). 

Below we present another type of nonstationary states, which are quite different from 
the above generalized coherent states. Recall that the problem was reduced to solving the 
equation ( 2.65 ) under the condition ( p.63|) . All the integrals of motion for such an equation 
can be constructed as functional combination of the operators 



fai, gaf 



02, 



''2 ' 



(2.87) 



whenever f,g obey the relations (p.69|) , ( p.70|) . Constructing integrals of motion that are 
linear combinations of these operators, we get coherent states. Any linear combinations of 
the operators ( |2.87|) do not commute with the operator (|2.8|) or Jz (|2.4|) . Thus, coher- 
ent states with definite values of these quantities cannot be constructed. The generalized 
squeezed coherent states ( p.44| ) and (|2.77| ) are eigenvectors of the operators A/i, A/2 (that fol- 
lows from ( p.79| ). The latter operators are integrals of motion and are quadratic in creation 
and annihilation operators. 



ATi = {At ^'^ - 4) - z,) , Af^ = {At "'^ - 4) (A^'^ - Z2 



(2.88) 



The operators Ai^ are defined in ( p. 68 ), and Ag''^ are defined in ( p. 35 ). The operators (|2.88 ) 
do not commute with Lz , Jz as well. 

One can see that besides the operators CL1CL2, at (^t 

, the only one quadratic combination 



that commutes with Lz , Jz is 



A = faia2 + gat^^t ■ 



(2.89) 



It is known |T5| that eigenvectors for such an operator can be normalized only for |/| > \g\, 
or for I /I = \g\. In the first case the eigenvectors have a finite norm, and in the second case 
they can be normalized to a 5— function. Let us consider the case |/| > \g\ only. Here the 
operator ( p.89| ) differs from 



AP 



aia2 



-2 + 4- 



p = pe — 1 < p < 1, K = uu^ + kq, kq = const (2.90) 



by a complex factor only. Thus, it is enough to consider the latter operator only. Let us 
demand that functions ip (u°, p, ip) be solutions of the equation ( |2.65| ), and, at the same time, 
eigenvectors of the operators A^, Lz, 



g,l 



I = 0,±1,±2, 



(2.91) 
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Such solutions can be constructed in terms of the Laguerre functions In,m (x) with noninteger 
indices, 



r,,i («°, P,^)=N exp (Uif - r) (1 + p)-" (1 - p)-^ J|,|+,,, (x) 
luu^ 1 + p"^ p — q 



= 1- . 
^ 2p ' 



|/| + 1 



X 



2pp 
\ — pi 



r(l+n) exp(-f 
^ r (1 + m) r (1 + n — m 



2p 



-X 2 $ (— m, n — m + 1: x) 



(2.92) 



Here $ (a,/5;x) is the degenerate hypergeometric function (see [|T4|, 9.210). For p^ = 1, the 
operator ( p.90| ) is anti-Hermitian and q is imaginary, Re g = 0. For p = 0, solutions have a 
very simple form 



u°,p,Lp] = Noexp{ilLp + q-To) Jm[2y/qp] , To = ^luJu" + ^, qe 



,0 I P ^^-ik 



(2.93) 

where Ju{x) is the Bessel function (see [|l4l, 8.402). The functions ( p.93| ) can be obtained 
from (|2.92|) as a limit p 0, as can be seen with the help of the property 



lim /, 



r+a,r+f3 



X 

4r 



Ja-I3 [x) 



(2.94) 



The functions (|2.92|) and ( p.93|) are orthogonal only with respect to quantum numbers /, 



Q 



',^„ijli)=ki'QF'{-s,-s'*;l + \l\-y), y=i^j 

r(i + |/| + s)r(i + |/| + s 



2p 



+ p' 



p^T{l + s)T{l + s'*) 



TtNN* l+Hi , 



(2.95) 



Here F(a,/3;7;x) is the hypergeometric function (see 9.100), and /^(x) is the Bessel 
function of imaginary argument (see [0, 8.404). Calculating (|2.95| ), we have used the 
integral table (see 0, 6.633.2; 7.622.1). 

The states ( |2.92| ) are not coherent states, however, they are, in a sense, close to such 
states. Indeed, let us consider the equations (pl6|) on classical trajectories. Then we get a 
classical relation 



P = P 



= Y L^h ^ + 4 |aia2|^ — aia2 — af 



(2.96) 



For p = 0, it follows from ( |2.91| ) that 0102 = —q, = hi. Thus, we can rewrite ( |2.96| ) in 
the form 



(2.97) 
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Calculating the mean value p by means of the functions ( p.93| ), we find 

P = Po + g + r, po = |/|-l-2|g| ^l'l~^^^'^P ■ (2.98) 

Thus, the time dependence of p is classical. The only constant which can differ from its 
classical value is po • 



III. EXACT SOLUTIONS OF RELATIVISTIC WAVE EQUATIONS IN 

LONGITUDINAL FIELDS 

A. Definition of fields 

Consider here longitudinal electromagnetic fields, which have the form 

E = nE, H = n/7. (3.1) 



Here n is a unit vector, = 1. Suppose n is directed along axis. Then, the fields ( p.l|) 
obey the free Maxwell equations whenever 

E = e(x^,x'^), H = H (x\x^) , 

where E (x°,x^) , if (x^,x^) are arbitrary functions of the indicated arguments. Thus, the 
fields under consideration can be represented by potentials of the form 

Ao = Ao x^) , Ai=Ai (x\x^) , A3 = A3 x^) ,^2 = ^2 (x\x^) , 

E = ^0^3 - 93A0, H = d2Ai - diA2 . (3.2) 



Thus, the operators ( p.6|) do not depend on the electric field (on Aq, A3). Therefore, imposing 
restrictions only on the magnetic field, we can maintain the relations (|2.6| - |2.1l|) . For a 



uniform magnetic field ( p.l| ), the commutation relations ( |2.7| ) are still valid and we use 
the semi-momentum representation, where these operators act on different variables (see 



Lorentz equations have the following form 

mx° + Ex^ = 0, mx^ + Ex^ = 0; 

mx^ + Hx^ = 0, mx^ + Hx^ = 0, x% = 1 , (3.3) 
which implies the following first integrals of motion, 

{{x,f + (X2)') = kl {{xof - (xsf) = + {hf , (3.4) 
where ki is an integration constant. 
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B. Klein-Gordon equation 

Consider the Klein-Gordon equation in the fields under consideration. Representing the 
wave function as 

^ = if , x^) ^ (^x\ x^) , (3.5) 

we find 

(p2 ^ p2 _ ^ (^x\x^) = 0, (V^ - -m^ - kf) ^ x^) = . (3.6) 

Using the variables x, y, 1], ^ defined in ( |2.5|) and ( p.l2|) , we can rewrite the equation for 
{x^, x^) in the following form 

{e-dl-k[')4^{x\x')=o, K' = ^. (3.7) 

The operator (^^'^ — (?| — /cf ^ does not depend on rj . Thus, it is convenient to go over to 
the semi-momentum representation, 

1 °° 

{x, y) = —j= / e^^y^i) (x, k) dk. 

— CO 

Substituting the integration over k by an integration over r], we obtain 

oo 

i;{x,y) = — / e-*^> (e, v) d^, ^ = V2x - r]. (3.8) 



TT 

— oo 



The function ip is an eigenvector of the operator A/" ( p. 141 ). The latter operator commutes 
with the operators from Eqs. (|3.6|) . In the semi- momentum representation we can write 

Af^Jn = n^n V) = Un (0 * (^) , (3-9) 

where $ (77) is an arbitrary function of rj. It follows from (|2.8| ) and the first Eq. (|3.6|) that 

Vl + Vl = 2Ti^^M + /i^7, kl = 2h^^n + . 

Solution of the last equation ( p.6|) can be found for fields that admit separation of the 
variables x^. For example, let us choose the potentials in the form: |e| Aq = A (x^) , A3 = 
0, |e| E = —d^A. In this case, stationary solutions of Eq. ( p.6|) read 

^ x^) = e-'^^o'V (a;^) , x" + RX = 0, R (a;^) = (ko + Af - - k^ . (3.10) 

Thus, the functions (3^) being written in the semi-momentum representation take the form 

v|>„ = e-'^o-'V (x^) (0 $ (r/) . (3.11) 

The equation ( p.lO| ) for x can be solved exactly for the following choices of the function 

A(x^): 



18 



A(x^^ = ax, A (x^) = a exp , A(^x^^ 



a 



x-^ 



atanh (l3x 



, A (x^) = a tan , A (x^) = [3 coth (/?x^) . 

The corresponding exact solutions are presented in [Q. Demanding that the Klein-Gordon 
function be an eigenvector of the operator (|2.16|) , we get an equation for the function 
$ (r/) from ( |3.11|) , 



a^a2$ {rf) = s$ {r]) . 



(3.12) 



Thus, $ (77) = f/s (r^) , and L^^ns = h{n — s) ■ Keeping this in mind and doing the 
integral (|3.8|), we obtain 



-1- r 



-ikox 



X + iy 



'2n \x-iy 



In,s (a;^ + 



where In,s are the Laguerre functions. Solutions of the equation ( |3.12| ) have been analyzed 
above. 



C. Dirac equation 

Let us present the Dirac wave function in the form 



M 



h {m + Vo + V3- ikia2) (Pi - 1V2) {m + Vo + Vs- iki02) 

{Vi + 1V2) [(m - Po - V^) as - hai] h [{m -Vq- P3) ^3 - hai] 



where p is an arbitrary spinor, which can be fixed by supplementary conditions. Then the 
following equations take place 



Vo - Vl - - k') if (x°, x^) =0, k' = kl + ieE, (3.13) 
(Pi + iV2) {x\ x^) = hk.^^i [x\x^) , (3.14) 
(7^1-2^2)^-1 {x\x^) = hkii>i (x\x'^) . (3.15) 

As a consequence of Eqs. (|3.14|) and ( p.l5|) , we get 

aiip_i = —i^Jnipi , alifji = i^/n■^|J_l, kl = 2'yn. 

Thus, we see that Vi = ipn-i , 'ip-i = —ii^n, and the problem is reduced to solving the 
equation (|3.13| ). The latter coincides with the second equation ( ^^ . 

Considering, for example, the constant and uniform magnetic field ( p.l| ) together with 
the electric field described by potentials |e| Aq = A (x^) , A^ = ^ \e \ E = —d^A, we get 

ip (x°, x^) = exp (-i/cox") X (x^) , x" + {ids + ko + A) x = ■ 
All possible solutions of the latter equation for the function x are presented in 
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PECULIARITIES OF INTEGRATION OF LINEAR DIFFERENTIAL 
EQUATIONS WITH NONCOMMUTATIVE SYMMETRIES 



Here we are going to return to the above results from a point of view of general theory of 
differential equation. Recall that we succeeded to find explicitly the transformation (|2.10|) - 
( p.l3| ) which has reduced effectively the number of the variables in the initial equations. In 
fact, that was the main starting point for all further constructions. However, one can see 
that this "reduction" of variables is a particular example of a general situation, which is 
described briefly below. 

Consider first the case of an integrable classical 2A^-dimensional Hamiltonian system 
with the Hamiltonian H . Suppose this system has independent integrals of motion that 
are in involution. It is well known that in such a case the variables of the type action-angle 
(J, are available, and the Hamiltonian depends on the action variables only, H = H{J). 
Let us suppose that for such a system exists one more independent integral of motion Y. 
Since Y is independent, it cannot commute with the former integrals, and, therefore, Y must 
depend on the angle variables. One can demonstrate that in such a case the Hamiltonian 
system is degenerate, det\\dH{J)/dJidJj\\ = 0, and, therefore, the Hamiltonian does not 
depend on some combinations of the action variables. For example, suppose the integral 
Y does not commute with the integral Jn only. Then the Hamiltonian can depend on 
the variables Ji, . . . , Jn-i only, otherwise H cannot commute with Y . Thus, we see that 
the noncommutative algebra of integrals of motion allows one to find canonical variables 
such that part of the corresponding action variables disappears from the Hamiltonian. This 
phenomenon is closely related to the topological properties of orbits for the Hamiltonian 
system. Namely, trajectories of the integrable Hamiltonian system with iV-dimensional 
commutative set of integrals of motion form (in the compact case) a winding of A^-torus 
in 2A^-dimensional phase space. If the set of the integrals is noncommutative then the 
dimension of the corresponding torus is r < (see ||16|| ). 

The phenomenon of variable "reduction" takes place in the quantum integrable Hamilto- 
nian system as well. As will be demonstrated below, by constructing a special isomorphism 
of linear functional spaces, we can transform the initial differential operator of an equation 
into another one with a reduced number of variables. The method which we are going to 
use for such a demonstration, is, in fact, the harmonic analysis on the noncommutative 
functional algebras. 

Suppose a differential equation 

H{x, d^)^{x) =0, tfjECc C°°(i?^) (3.16) 

with N independent variables x G i?^ admits a noncommutative algebra of functionally in- 
dependent symmetry operators JF = {Xa{x, d^)}- The corresponding commutation relations 
are in the general case nonlinear, 

^[Xa,Xh] = nab{X), a,b =l,...,n = dim J^.. (3.17) 

Here Qabi^) are symmetric operator functions. The linear case, when Qabi^) = C^^Xc, 
corresponds to a Lie algebra, the quadratic symmetric functions Qab{X) correspond to a 
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quadratic algebra, and so on. The algebra T corresponds to the algebra T' = {Ya{x,dx)} 
of the invariant operators on C: 

[Xa,Y^]=0, ^[Y^,Yfs]=iu^piY), a,P = l,...,n' = dmir. (3.18) 

We denote via E{J^) and E{J^') enveloping fields for the algebras and T' respectively. 
Elements of E{J^) and E{T') are symmetrized operator functions of the generating operators 
Xa Ya- It is clear that the centers of the enveloping fields coincide, i.e. Z{E{J^)) = Z{E{T')). 
The elements of the center Z = Z{E{J^)) are called Casimir operators. The number of the 
independent Casimir operators, which generate the center Z is called the index of the algebra 
T: r = indJF = indJF'. If we replace the operators X and Y in the operator functions 
1^0/3 (^) and Q.ah{X^ by arbitrary complex numbers ^ and /, then the index of the algebras 
JF and JF' can be calculated according to the formula 

r = sup coranknafe((^) = sup corankco'Q,^(/) . (3.19) 
Cec /ec 

One can show that the following relation takes place 

n + n' = 2iV. (3.20) 
Let us introduce the notion of the A-representation of the algebra |r7|. In fact, the 



A-representation is the result of the quantization of the classical Poisson bracket and can be 
understood as a realization of the algebra T by an irreducible set of differential operators 
X = X{q,dg,j), dependent on r parameters j = {ji, . . . , jr), and acting in a space of 
functions of [q] = {n — r)/2 independent variable^ q & Q , 



^[Xa,X,] = -QabiX), K^{X{q,d„j)) = ^^{j), det 



^0. (3.21) 



Here iiT^ are all the independent Casimir operators of T. In a similar manner, we construct 
the A-representation {Y} of jF'in a space of functions of [q'] = {n' — r)/2 independent 
variables q' E Q' , 

'-[Y^.Yp] = cu^^iY), K'^iY{q',d,,,j)) = K,iX{q,d„j)) = k,{j) . (3.22) 
Suppose that in the spaces of functions of x, of q, and of q' are defined scalar products 



(<^,V')= J <^{x)ij{x)dfi{x), {ip,%l)) = j (f{q)ij{q)dfi{q), {^,^y = J ip{q')ip{q )dii{q), 

RN Q Q> 

(3.23) 

where d^{x), dfi{q), and d^{q') are some measures on i?^, Q, and Q' respectively. And 
suppose that the operators Xa{x,dx), Ya{x,dx) and the operators Xa{Q,dq,j),Y{q',dq',j) 



^Via [q] we denote the number of the variables q. Similar notations are used below. 
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are self-conjugate with respect to the corresponding scalar products^. Now we define the set 
of distributions Di,( solution of the overdetermined system of the equations: 

Xa{x,d^) - Xa{q,dg,j) 

The distributions obey the completeness and orthogonality relations: 

J Dlg,{x)Dlg,{x)dfi{x) = 6{j,j)5{q,q)5{q',q); (3.25) 



X 



Yaix,d^) -Ya{q',dq>,j] 



^qq' 



0. 



(3.24) 



Diq>{,x)Dl ,{x) dn{i)dn{q)dn{q') = 5{x,x). 



(3.26) 



Here dfi{j) is the spectral measure of the Casimir operators K{X) (= K'{Y)). Due to Eqs. 
( p.21| ), (|3.22|) ( 3.24 ) the distributions Dgg,{x) are eigenfunctions of all the Casimir operators, 



K^{X{x,d^))DL{x) = K^{j)Dlg,{x), /i = l,...,r. 



(3.27) 



Usually one can find D-j^^, (x) by an integration, at least in the case when ^ is a Lie algebra. 
As a consequence of ( p.26[ ) and ( p.25[ ) we can define the direct and the inverse Fourier 
transforms: 



^(g,g',j) = J Dlg'{x)^{x) dn{x) , 

^ix)= / Dlg,{x)ip{q,q'J)dfx{j)dfx{q)dfx{q'). 



(3.28) 
(3.29) 



The Eqs. ( p.28| ) and ( p.29| ) establish an isomorphism of the spaces C and C = {ip}. It is 
important to stress that under such an isomorphism the operators X, Y are transformed 
into the differential operators X, Y that act in the spaces of functions, which depend on a 
fewer number of variables, 



X{x,d^)^/j{x) ^ X{q,dg,j)ilj{q,q',j), Y{x, d^)ilj{x) ^ Y{q',dg>J)ip{q,q',j) 



(3.30) 



Let us return to the equation ( p.l6 ). Here we can conclude that H G E{J-'') since the 
operator H commutes with all the operators of the algebra J-' = {Xa}. In turn, that means 
that there exists an operator function H{Y) such that H{x, d^) = H{Y{x, dx))- Let us look 
for solutions of the equation (|3.16| ) in the form (p.29|) . Using the isomorphism (|3.30| ), we 
get: 



H{Y{q',dg,,m{q,q',j] 



0. 



(3.31) 



Thus, departing from the equation ( p.l6| ) without loss of any information, we have arrived 
to a differential equation with N = [q'] = {n' — r)/2 independent variables. Taking into 
account ( ^.20] ), one can obtain 



This supposition is not necessary and is introduced to simplify the consideration. 
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iV = iv_ l(dimJ^ + indJ^). (3.32) 

Thus, the existence of a noncommutative symmetry algebra results into the phenomenon 
of variable reduction. Indeed, we started with N = [q'] + [q]+r variables, [q] variables have 
disappeared completely, and r = [j] variables remain in the equations as some parameters. 
Respectively, the solution of the equation ( |3.31| ) contains as a factor an arbitrary function of 
the variables q,j. The number of variables q is equal to [g] = (n — r)/2 = (dim JF — ind jF)/2. 
In the commutative case: indJF = dimjF, and [q] = 0. In the noncommutative case: 
dimjF > indJF, and [q] > 0. Thus, the reduction of variables always takes place when there 
exists a noncommutative algebra of the integrals of motion. 

Let us apply the above consideration to the Klein-Gordon equation ( p. 21 ) in a uniform 
magnetic field. In this case we have four {N = 4) variables and five (n = 5) independent 
symmetry operators: JF = {Pq, -P3, 0,2, at^ -^2 = fi-L}, 

Pq = ihdo, P3 = ihd^, L = udu — udu , a2 = du + u/2, = —du + u/2 , u = x + iy . 

The non-zero commutation relations are [02,02'] = 1; [-^) '^2] = 0-2, [L,a2] = —at- It 
follows from ( |3.19| ) that r = indJF = 3. = [q'] = JC E Z according to ( p.32| ). Thus, the 



equations ( p.31|) present algebraic relations on the parameters j (and on the parameters of the 
equation itself ). Besides, dim^' = 3 = r, due to ( p.2(]| ). Thus, the algebra of the invariant 
operators is placed completely in the center. Or more simply, there are no operators Y, Y 
and variables q' in the case under consideration. The center Z is generated by three Casimir 
operators, those are Ki = Pq, K2 = P3, K3 = J\f = L + a^a2 = L + |(a^a2 + 020^ — 1). 
Let us construct the A-representation of the algebra JF: 

^0 = ji = , P3 = j2 = ^A^s , 5.2 = q, at = dg, L = -qdg + n, n = = 0, 1, . . . . 



The operators 0.2 and at are mutually conjugate with respect to the scalar product (p.23|) 



with the measure dfj,{q) = exp{—qq)(Pq/7i, {(Pq = dqidq2, q = qi + iq2)- The operator L is 
self-conjugate, and the space C is built up from analytic functions dependent on the variable 
q and on the parameters j = {kQ,k3,n). Here the Casimir operator A/" has the following 
form Af = L + ^{dtd2 + 02^2^ — 1) = n. 

From the equations p.24| ) we can find the set D^g{x) that obeys the completeness and 



orthogonality relations. Such a set has the form 

P)i(a;) = e*(fcoxO+fc3x3)ggs-u«/2^^ _ g)"/(27rv^) , (3.33) 
J Dl{x)Dl{x)dx = 6iko - ko)5iks - h)mQ)Snn, (3.34) 



n=0 



J2 / Di{x)Di{x) dkodk^dniq) = 6{x^ - x°)(5(x^ - x^)6{x - x)6{y - y). (3.35) 



In Eq. ( p.34| ) S{q,q) = exp(gg) is a ^-function with respect to the measure dfi{q). To justify 
the validity of (|3.34| ) and ( p.35| ) one can use the following relations |T^ 

00 

Vn{q)vm{q) dn{q) = 6nm, =^{q,(1) =exp(gg), Vn{q) ^q^-fVnJ.. 

n=0 
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As was already mentioned above, the Klein-Gordon operator K, belongs to the center Z 
, and , therefore, can be presented in a polynomial form of the Casimir operators Pq, PsyN", 
which generate this center. Such a representation is given by Eq. ( |2.9| ). The Klein-Gordon 
equation in the space C, i.e. the equation (|3.31 
parameters j = {ko,k3,n). Thus, the functions (| 
equation (12 



is, in fact, the relation ( p.21| ) for the 
33f) form a basis of the Klein-Gordon 



(with allowance made for (2.21 
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